5*0 pairing correlation in symmetric nuclear matter with Debye screening effects 
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The ^5*0 pairing of symmetric nuclear matter is discussed in the frame work of relativistic nuclear 
theory with Dyson-Schwinger equations(DSEs). The in- medium nucleon and meson propagators are 
treated in a more self-consistent way through meson polarizations. The screening effects on mesons 
due to in-medium nucleon excitation are found to reduce the ^Sq pairing gap and shift remarkably 
the gap peak to low density region. 
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Compared with nonrelativistic frame work, the rela- 
tivistic nuclear theory can successfully describe the sat- 
uration at normal nuclear density. The basic meson 
exchange is normally considered as the nuclear satura- 
tion mechanism. The original a — lo theory of quantum 
hadrodynamics model (QHD) developed by Walecka et 
al. and its various extensions have been widely used 
to discuss the properties of finite nuclei and nuclear 
matterHHUll. 

Superfluidity of strongly interacting Fermi system is 
very important for understanding the properties of finite 
nuclei, such as the dramatic reduction of the moments of 
inertia in rotating nuclei or the energy gap in the spectra 
of many even-even nuclei HQ. The existence of superflu- 
idity may also affect the dynamical and thermal evolution 
mechanism of neutron stars because it is closely related 
to the emission of neutrino and cooling of neutron-rich 
matter. It is also argued that the superfluidity of nuclear 
matter can lead to the glitches of astronomy phenomena 
and attracts much attention in contemporary physics 7] • 

Although there are many works in the literature ad- 
dressing the superfluidity of nuclear matter, the main re- 
sults are obtained from the non-relativistic nuclear the- 
ory and no definite conclusion can be made yet. We 
noted that since K. Kucharek and P. RingQ first de- 
rived the relativistic Hartree-Fock-Bogoliubov equation 
by using Green function method and the Gor'kov fac- 
torization analogously to nonrelativistic BCS theory, it 
was found that the superfluidity gap value is about three 
times larger than the "standard" value obtained with the 
nonrelativistic Gogny force To improve the descrip- 
tion of superfluidity with relativistic nuclear theory, var- 
ious approaches have been investigated, such as using 
external potential as input or various cut-offs of the inte- 
gration momentum [iCl. UlL . To our knowledge, there 
is even no definite result about the superfluidity gap(the 
gap value and peak position) for symmetric nuclear mat- 
ter within nonrelativistic or relativistic nuclear theory up 
to now. However, it is widely accepted that the ^5*0 gap 
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values at normal nuclear density should be very small. 

Essentially one can not expect that the softness of 
equation of state(EOS) describing the bulk property of 
nuclear matter is directly related with the superfluidity 
property of nuclear matter. For example, although the 
nonlinear a-Lu model with the possible embarrassing neg- 
ative coupling constants b and c (which in principle lead 
to instability of nuclear system at high density scenario) 
in a self-interaction terms ba^ + ca'^ can give a very soft 
EOS with mean fleld theory(MFT)[3, the gap behavior 
is similar to that in the original version by using frozen 
meson propagators 

Theoretically, since it is difficult to make low en- 
ergy calculations directly with quantum chromodynam- 
ics(QCD), one has to work with effective theories. As an 
effective theory, QHD-I model (and its extensions) has 
been widely used to discuss the effective meson masses 
under extreme environment in the past[ll[ll|il[l3. 
In principle, when one discusses the in-medium proper- 
ties of nucleons and mesons, one has to take into ac- 
count the back-interactions of nucleons with in-medium 
mesons. Therefore the resummed nucleon and meson 
propagators would form a closed set of coupled equa- 
tions and should be solved simultaneously. With this self- 
consistent way, a softer EOS with an acceptable compres- 
sion modulus K, in dealing with realistic nuclear matter 
can be obtainedpHTTj. In the spirit of mean field theory, 
the exchanged mesons in determining nucleon propagator 
are not free but medium dependent. Their masses should 
be determined together with the nucleon mass through 
Dyson-Schwinger equations self-consistently, as indicated 
by FigCl 

It would be very interesting to analyze the in-medium 
effect of mesons on the superfluidity of nuclear matter 
in the frame work of relativistic nuclear theory. In su- 
perfluidity state and with QHD-like Lagrangian, Dyson- 
Schwinger equations for the nucleon and meson propa- 
gators as indicated in Fig^and the energy gap equation 
as indicated in Fig|21 form a new closed set of coupled 
equations. The in-medium meson propagators D instead 
of the normally used bare ones will affect the ker- 
nel in the BCS gap equation. We will see below that 
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FIG. 1: Diagrammatic representations for the propagators of 
in-medium nucleons(a) and mesons(b). 

the gap behavior with in-medium meson propagators is 
quite different from that with bare ones. We found that 
the polarization effects leading to screening have been 
widely discussed in the nonrelativistic frame work of nu- 
clear theory [11 H HE IH H^, [13 , this problem has not 
yet been addressed before in terms of relativistic field 
theory. In this letter, we want to discuss the effects of 
in-medium effective potential for nucleon-nucleon inter- 
action on the superfluidity due to screening by using the 
original renormalizable formalism of cr — a; model. 

Let us start with the four-dimensional gap equation by 
using the standard Nambu-Gor'kov formalism in the lad- 
der approximation of the meson exchanges as indicated 

by FigEiillillii, 

^ < ^'^'^ > ^^^^^ 

where K is the four momentum K = (fco, k), < P|r|_R' > 
is the interaction kernel and 
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is the Nambu-Gor'kov anomalous propagator with 
s{K) = Ek — Ekf being the quasi-particle energy above 
Fermi-surface. For ^5*0 pairing, the gap equation can be 
reduced to|i.lTot 

(2) 

with Ek ^ E* + X and E* = y^M*f+^. The quan- 
tity A related with the baryon current is obtained from 
the tadpole self-energy of nuclcon propagator with in- 
medium vector meson in Fig^ 

2 roo 

^ = ^T2 ^l^^dK (3) 



where 7 = 4 is the spin-isospin degeneracy factor for 
symmetric nuclear matter, and vf, is the BCS occupation 
number 




FIG. 2; Diagrammatic representations for the gap equation 
and the interaction kernel in instantaneous approximation 
with screened meson propagator. 

The interaction kernel < PlFjif > in our treatment is 
approximated by the in-medium meson propagators in- 
stead of the bare ones. We use the static(instantaneous) 
approximation by neglecting the retarding effects 
Since the meson propagators Z?CT.t»j(0, k) with the van- 
ishing temporal component of four-momentum are now 
medium dependent, the Debye screening effects will play 
an important role in the in-medium particle-particle in- 
teraction potential 

f;(p,k) =< ps', ps'|y|ks, ks > - < ps', ps'|y |ks, ks > 
, ^ M*^ Tr[A+(k)rA+(p)7°T+r+T7°] 
^2E*{k)E*{p) + '^^ 

where A+(k) = ^2M*" projection operator of the 

positive energy solution and T = 47^7'^ is the time re- 
versal operator. The tilde typescript "~ " means the 
Gor'kov time reversal state and F is the corresponding in- 
teraction vertex oi a /lo with nucleons. The assymetrized 
matrix elements Vpp{p,k) in the gap equation Eq.Q for 
^5*0 pairing is obtained through the integration of ^(p, k) 
over the angle between the three-momentums k and p 

Vpp{p,k) — J w(p, k)dcos6'. (6) 

The effective nucleon mass is determined by the rel- 
evant mass gap equation through tadpole self-energy of 
nucleon propagator with scalar meson 
2 p 00 /\ //" * 

mi 27r^ Jo E; ^ 
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with AMyl^ being the vacuum fluctuation contribution 
9l 1 
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The polarization tensors Il^^^{kQ,]i) determining the 
in-medium a and uj propagators are calculated by using 
corresponding Dyson-Schwinger equations as shown in 
FigH For brevity, here we list only the sigma meson 
self-energy explicitly, 
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with 



= In 



fc2 - 2p|k| - 2koE; 
fc2 + 2p|k| - 2kQE* 



The effective masses fha and ifi^ in Eq.Q and Eq.jTI) 
are determined by the corresponding polarization tensors 
with vanishing four-momentum transfer, and the Debye 
screening masses to* and to* in assymetrized matrix el- 
ements Vpp{p, k) are determined by the pole positions of 
corresponding spacelike propagators ^(0, k) due to 
taking the static approximation 16J. For example, the 
transverse mode screening mass to* is determined self- 
consistently by 



TO:^ = TOi + n^(o,zTO:), 



(10) 



where is the transverse part of polarization tensor 
n^''(fco, k). In principle, the longitudinal mode screening 
mass is different from the transverse mode one for in- 
medium vector meson due to the broken Lorentz invari- 
ance. However, neglecting this little difference doesn't 
affect the qualitative result in realistic numerical calcu- 
lation. 

Considering the in-medium meson effects on the prop- 
erty of nuclear matter, one should refix the parameters in 
the model. Noting that the effect of superfluidity on the 
bulk property is negligible, we fix the parameters by nor- 
mal nuclear matter with saturation condition of binding 

energy e„ = —15.75 MeV at the normal nuclear density 
lo 



with /cf = 1.42/to ^. The relevant parameters are listed 



in TablelD 

The remaining task will be the numerical solution of 
the coupled equations indicated by FigsQ] and [3 It 



should be noted that the relativistic kinematic factors 
guarantee the convergence of the gap equations such as 
Eq.(j2l for the relativistic nuclear theory and lead to a 
definite result for the gap. In principle, the momentum 
integration upper bound in relevant equations such as in 
Eq.(|2Jl is infinity. However, a concrete upper bound must 
be used to give a numerical result by solving the gap in- 
tegral equation. Strictly speaking, the gap value should 
not be sensitive to the adopted momentum upper bound 
and the sensitivity of momentum cut-off on the gap has 
been analyzed in such as in Refs.|^ which can be 
also reflected by the gap function indicated by Fig. |2fb) . 
A concrete and large enough momentum upper bound 
Ap = 20/to~^ has been used in this work for the descrip- 
tion of screening effects. To focus on the characteristic 
due to polarization effects, the a — u mixing effects have 
been neglected, which will not affect the result qualita- 
tively although it deserves further study. 
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FIG. 3: (a) The pairing gap A/ at the Fermi surface as a 
function of density characterized by the Fermi momentum kf. 
(b) The gap function A(fc) as a function of momentum k for 
fixed Fermi momentum kf = 0.5/m~^. The dot-dashed lines 
correspond to the result obtained by MFT and dashed lines 
correspond to RHAT^, while the solid lines correspond to 
our self-consistent resummation approach. 



The numerical results of the superfluidity gap equa- 
tions are shown in Fig|3 In the upper panel (a), we 
indicate the gap curves A(fcj) versus Fermi momentum 
kj, and the gap functions A(/c) at given Fermi surface 
momentum kf = 0.5fm~^ are shown in Fig. Olb). 
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FIG. 4: Assymetrized matrix elements Vpp{k, kf) in momen- 
tum space at the Fermi momentum kf — 0.8/m~^. The hne- 
styles are similar to those in Fig0 

Compared with the previous superfluidity results of 
relativistic nuclear theory in the literature, the gap value 
we obtained is very small and the peak position is shifted 
to the low density region remarkably. As mentioned in 
the introduction, this interesting result is not due to the 
softener EOS but much attributed to the screened ef- 
fective particle-particle interaction potential as indicated 
by Fig21 The key point is that the a and cu propagators 
in the gap equations are not bare but in-medium ones 
determined self-consistently by Dyson-Schwinger equa- 
tions. The effective nucleon-nucleon interaction poten- 
tial with Debye screening of in-medium mesons leads to 
the change in the interaction force range. Different from 
the scenario of bare meson propagators used in the gap 
equation, the particle-particle potential is more sensitive 
to density, which can be understood from the correspond- 
ing attractive and repulsive force range changes for differ- 
ent densities characterized by Fermi momentum fc^(not 
displayed obviously in Fig. 

It is clearly demonstrated in Fig|31 that our self- 
consistent approach reduces the difference between the 
nonrelativistic and relativistic theories about the max- 
imum gap value and the peak position of ^Sq pairing 
correlation. The significant improvement by the self- 
consistent resummation approach for the particle-particle 
interaction leading to pairing is reflected on two aspects: 
One is at the saturation density with fc^ = 1.42/m~^, 
the other is at kf = 0. The improvement at A: = is 
crucial by noting that the MFT and RHA approaches 
with frozen meson propagators give unrealistic non-zero 
gap values - 1.94 MeV/0.36 MeV[l3. 

Summarizing, with a set of more self-consistent equa- 
tions for the resummed in-medium nucleon and me- 
son propagators and the superfluidity gap by Dyson- 
Schwinger Green function approach, we have studied the 
^5*0 pairing correlation in symmetric nuclear matter and 
compared our results with those obtained by MFT and 
RHA approaches. The Debye screening effects of in- 
medium meson propagators can reduce significantly the 
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TABLE I: The parameters determined at normal nuclear 
density with k°f = 1.42/m"^ in MFT, RHA and our self- 
consistent approach(labeled as SA) with Mat — 939 MeV, 
niu, — 783 MeV and = 550 MeV. We show also the com- 
pression modulus K, (in MeV), the medium dependent coupling 

constants (determining EOS) ^ and Cl = gl 

the maximum of gap value A™ (MeV), the peak position fc™ 

{fm~^) and the "gap value" A(0) (MeV) at kf = 0. 



2 2 

ga g^ 






fC 




Af 


kf 


A(0) 


MFT 91.64 136.2 


267.11 


195.87 


545.43 


0.556 


9.4 


1.0 


1.94 


RHA 62.89 79.78 


183.31 


114.73 


468.24 


0.718 


8.3 


1.0 


0.36 


SA 48.90 53.40 


123.17 


66.078 


338.00 


0.794 


1.9 


0.5 






superfluidity gap value, while the gap peak position is 
shifted remarkably to low density region. 
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